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Abstract 

In this paper we study the quantum cosmological Kantowski-Sachs model and solve the 
Wheeler-DeWitt equation in minisuperspace to obtain the wave function of the correspond¬ 
ing universe. The perfect fluid is described by the Schutz’s canonical formalism, which allows 
to attribute dynamical degrees of freedom to matter. The time is introduced phenomenolog¬ 
ically using the fluid’s degrees of freedom. In particular, we adopt a stiff matter fluid. The 
viability of this model is analyzed and discussed. 
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1 Introduction 

At first sight it may seem meaningless to study quantum cosmological models msiEiiii taking 
into account only an absolutely finite number of degrees of freedom, i.e., a minisuperspace. In 
general, the idea is to restrict the number of degrees of freedom of the metric to a finite number, 
by imposing determined symmetries. As a consequence we have that the quantum constraints 
become easier to handle, when written in the usual variables; since we allow fewer metric varia¬ 
tions when getting the equations of motion from the action, we may get more classical solutions, 
some spurious ones being unstable; from both semiclassical and quantum point of view, the 
results may not be equivalent to appropriate sectors of the full theory. Despite this, the con¬ 
struction of a quantum theory of gravity must be related, among other things, to obtaining the 
solution to the Wheeler-DeWit equation in the superspace. Hence, to exclude infinite degrees 
of freedom does not seem to be a reasonable procedure. However, the quantization in minisu¬ 
perspace is still very useful for other reasons. Minisuperspace quantum models are very simple 
systems, without much of symmetry of general relativity. They provide non-trivial examples 
where we can apply the ideas of quantum gravity and mathematical techniques. They are still 
relevant for the discussion of major conceptual issues of quantum gravity such as the problem 
of time and the invariance under time reparametrization of classical general relativity and the 

*e-mail: f lavio. alvarengaOuf es .br 
A-mail: rfracalossi@gmail.com 
G-mail: rodolfo . camargoSpq. cnpq.br 
R-mail: sergio. vitorinoSpq. cnpq.br 


1 



physical interpretation of the Hamiltonian constraint. In some cases, a time coordinate may be 
identified with the space volume, which is a growing function in an expanding universe. But, all 
these attempts have limited applications until now, and the problem of time in quantum gravity 
and quantum cosmology remains an unsolved puzzle la 0 El El E]. Even sometimes working 
with only the existence of the gravitational field, quantum aspects can be consistently verified 
or at least induced. This is even more evident in the study of quantum properties of black holes 
as well as the use of the known quantization processes to the universe as a whole. 

Focusing only in minisuperspace we can ask what are the degrees of freedom responsible for 
ensuring certain fundamental aspects of the desired description of the cosmological model. Only 
after verifying these aspects we can take into account the thermodynamic aspects of the problem 
and the presence of other fields. However, the use of metrics with additional degrees of freedom 
to those found in FRW metric can be deeply complex, especially from a mathematical point of 
view. As a result, models that represent inhomogeneous and anisotropic universes are largely 
unexplored in comparison to the FRW. We can cite Bianchi models, for example jlUj . These 
cosmological models have a limited number of analytical solutions, and to obtain numerical 
solutions for the Wheeler-DeWitt equation can be a complex process. 

In this work, we study the quantization of minisuperspace models within the context of the 
Kantowski-Sachs metric, and we describe matter as a perfect fluid with an arbitrary barotropic 
equation of state p = ap. Therefore, we will adopt Schutz’s canonical formalism mm, which 
describes a relativistic fluid interacting with the gravitational field. This description is essentially 
semiclassical m, but it has the advantage of assigning dynamical degrees of freedom connected 
with the matter, which can naturally be identified with time, leading to a well-defined Hilbert 
space structure. In particular, the stiff matter fluid is used. The importance of the stiff matter 
in cosmological models is discussed in m- 

The Kantowski-Sachs space-time m is a cosmological model technically less complex than 
other anisotropic metrics. It represents a homogeneous and anisotropic geometry. One may 
question the relevance of anisotropic models of the observational point of view since the universe 
today seems to be homogeneous and isotropic. However, there is no reason, a priori, that the 
universe has always been so. In early times, isotropization processes of quantum nature or even 
thermodynamics nature m must have occurred on an essentially less regular space-time. In 
this sense, the study of anisotropic models in the early universe where quantum effects have 
great relevance is justified. 

The Kantowski-Sachs space-time also has some interesting features. First, its classical so¬ 
lutions are known in different contexts, as well as their quantum solutions da EZl Eg [H]. 
Furthermore, Kantowski-Sachs metric may coincide with the vacuum Schwarzschild solution 
inside the singularity. As it may also describe the interior of black holes, it can also be used 
to discuss their quantum states and its entropy. Models such as the Kantowski-Sachs are still 
interesting in order to study the behavior of the added degrees of freedom in quantum cos¬ 
mological models. Finally, it is possible that the construction of a Kantowski-Sachs quantum 
cosmological model suggest adaptations and modifications in the quantization methods applied 
to cosmology. For example, in j2D], the quantization in minisuperspace of the Bianchi I model 
reveals a startling discrepancy of interpretations of Bohm-de Broglie and many worlds of quan¬ 
tum mechanics: the quantum model is non-unitary and the norm of the wave function becomes 
dependent of time. Perhaps the analysis of other anisotropic quantum models is likely to explain 
this anomaly. In general, the condition that an operator is self-adjoint is sufficient to provide a 
spectral behavior, as unitarity, appropriate to a quantum system. The appearance of anomalies 
in functions associated with a self-adjoint operator must indicate that there are other relevant 
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mathematical aspects that must be taken into account in the construction of a quantum cosmo¬ 
logical model and that go beyond spectral analysis. In particular, there is suspicion that even 
the build procedure of the wave packet may introduce such anomalies. 

Here we apply the previous investigations of quantum perfect fluid models isniEHEaisiiMi 
ES] to quantize the Kantowski-Sachs model. The Wheeler-DeWitt equation is solved for the min¬ 
isuperspace and the wave function is obtained. In this procedure the time is phenomenologically 
introduced through one of the degrees of freedom of a stiff matter fluid. 

This article is organized as follows: In the next section, we show the Kantowski-Sachs metric 
and calculate the Hamiltonian function from the gravitational part. Then, the quantum cosmo¬ 
logical model is constructed using the approach of the phenomenological time. In section we 
proceed a detailed discussion and conclusion of the results obtained. 


2 The cosmological model: the choice of time and the Wheeler- 
DeWitt quantization 

In this model we consider a flat, homogeneous and anisotropic universe described by the Kantowski- 
Sachs metric. According to m the Kantowski-Sachs spacetime is not a spacetime itself, but 
only part of it. This restriction holds for the vacuum model, but disappears if a perfect fluid 
is coupled to gravity. Here, we introduce a perfect fluid, described by the Schutz’s canonical 
formalism mm to obtain an explicit time variable related with the fluid degrees of freedom, 
solving two problems: the first one, described in j2B], and the problem of the time variable that 
is absent in the quantum cosmological models. 


2.1 Kantowski-Sachs spacetime 


We start with the Einstein-Hilbert action plus a boundary term 

Sg = [ d'^Xy/^R + 2 [ d^xVhK , (1) 

JM JdM 

where R is the scalar curvature, taken as a function of and its derivatives, K = habK°'^ such 
that Kab is the extrinsic curvature and hab is the induced metric over the three-dimensional 
spatial hypersurface, which is the boundary dM of the four dimensional manifold M. The 
factor IGttG is made equal to one in the first term of equation Q, where G is the Newtonian 
gravitational constant. 

We will consider a model of a universe that is homogeneous but anisotropic. Thus, we shall 
deal with a spacetime such that the line element can be put in the more usual form m 

ds^ = —N{t)‘^dt^ + a{t)‘^dr‘^ + b{t)‘^{d6‘^ + sin^ 6 dip^) , (2) 


where a{t) and b{t) are the scale factors and N{t) denote the lapse function. This is known as 
Kantowski-Sachs metric m- 

We can compute the Ricci scalar R, which for the Kantowski-Sachs metric turns out to be 

a A h N . N b 4 db 2 6^ 2 

l^ab^ 

The extrinsic curvature K is given by 


K = - 


1 

N 



(4) 
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Replacing the equations Q and @ in Einstein-Hilbert action 0 we have the following 
expression 

^7 /r,A 7 4a66 2ab‘^\ 

dt I 2Na -- 7 ^ 1 . (5) 


83 = 


N 


N 


There are two types of singularities in this geometry: cigarlike singularities, when 6 —)• 0 and 
disklike singularities, when a —)■ 0 m- However, the curvature of the hypersurface = 
divergent for the cigarlike singularity. Due to the absence of matter, the disklike singularity is a 
singularity where the right handed coordinate system becomes a left handed system. Moreover, 
it is an indication for the incompleteness of the vacuum Kantowski-Sachs spacetime mm- 
The model may be completed by adding matter, e.g., in the form of stiff matter, like in this 
work. 

With the gravitational action (|^, the gravitational Lagrangian of the Kantowski-Sachs model 
is given as follows 

abb aiP Na 

^ ~ ~lv “ ^ ^ 

In order to deal with a mathematically tractable problem let us perform the reparametriza- 
tion 


a 


^ c ac 


Thus, the purely gravitational Lagrangian is written as 


^ Na ^ 

2Na^ 


c" 

2Na 


( 6 ) 


Using the canonical formalism the gravitational Hamiltonian can be expressed as 


IT “ r 2 2 2 2 2i 

Hg = i^[a p^-cp^-c]. 


(7) 


In the next section we will introduce the matter as a barotropic perfect fluid described by 
the thermodynamic variables. 


2.2 Phenomenological time: Schutz’s canonical formalism 

Schutz mm introduced, in the context of general relativity, a representation for a perfect 
fluid in terms of velocity potentials. The formalism describes the dynamics of a relativistic fluid 
in interaction with the gravitational field. 

The hydrodynamic equations are put in Eulerian form, with the four-velocity expressed in 
terms of five potentials /3,6 and S: 

Uu = + CP,u + ) ( 8 ) 

d- 

each satisfying its own equation of motion. In this expression p, is the specific enthalpy, the 
variable S is the specific entropy, while the potentials ^ and /3 are connected with rotation and 
are absent in this model. The variables (p and 9 have no clear physical meaning. 

The derived dynamical equations provide a description of the usual equivalent hydrody¬ 
namic equations based on the divergence of the energy-momentum tensor. Erom this, follows 
an especially simple action. 
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( 9 ) 


Sf = J (fxyf^p , 

where p is the pressure, which is associated to the energy density by the equation of state p = ap. 

It is important to mention that the Einstein’s equations for a perfect fluid can be found 
when we vary the action (© plus the Elnstein-Hllbert action In respect to the metric tensor. On 
the another hand, the changes in relation to the velocity potentials is similar to the Eulerlan 
equations of evolution of the fluid, 

= 17"C. = 0; [/"A. = 0; 


= = 0; (poC/");. = 0. 

Here Tg is the temperature and po is the rest mass density of the fluid. 

The basic thermodynamic relations take the form 

p = po(l+n); /U = (l + n) + —; TgdS = dU + pd(, 

Po \PoJ 

n being the specific internal energy. 

If we write 


dU 


+ pd (—) =(1+ Il)d[ln{l + n) — a Inpo] , 
\PoJ 


we can Identify = 1 + H and S = /n(l + H)/po“- 

After a few mathematical steps, the energy density and pressure reduces to 


( 10 ) 


( 11 ) 


( 12 ) 


P = 


P 


i+- 


^ H" 1 

By means of the normalization condition 


e “ and p = a 




Q: H” 1 


l+- 


oi _ 5 

e 


(13) 


U’^Uu = -1 , 


once can express p in terms of the potentials 


(14) 


P = ^{^ + SS) . 


(15) 


Elnally, by using (13) an (151 in the matter action ([9|), it is possible to obtain the Lagrangian 
of the fluid: 


Lf = — 


a 




(16) 


a (a + 1)1+1/° 

The fluid conjugated momenta are derived from the above Lagrangian, written in terms of 
the canonical variables 


P<f> = 

PS = 
Pe = 


— _ 

a (a+ 1)1/° 




0 


5 


(17) 








and again by the canonical formalism we obtain 


a 




This can be put in a more suggestive form by means of the canonical transformations 
T = , pr = 0 = (/> - (a + 1)^ , = P 4 ,, 


and the super-Hamiltonian (18) can be written in the final form: 


PT 


( 18 ) 


(19) 


( 20 ) 


Time introduced via Schutz’s formalism is a global phase time. In fact, to be a global phase 
time the function must be monotonously increasing, defined in phase space and satisfying the 
relation 

[ t , H ']>0 . ( 21 ) 


In particular, the time introduced in equations (^19^ and (20) is an extrinsic time. If the function 


defined in phase space is of type t{qi,pi), it is called extrinsic time. But, if the function has the 
form t{qi) it is called intrinsic time |5j. The idea of a global phase time in anisotropic universes 
of this type is well exposed in reference [IB]. It is possible, in some contexts, to connect the time 
with the geometric origin of the quantum theory. In the Kantowski-Sachs universe the time can 
be a global phase factor if, throughout some gauge transformations, we can make it obeys the 
relation [t,H] > 0. 


2.3 Kantowski-Sachs quantum cosmology with perfect fluid 

In this section, we describe the Schutz’s formalism in Kantowski-Sachs spacetime. The total 
Hamiltonian H = Hg + Hj is given by the gravitational part Hg plus the fluid part Hf, that is 



22 22 

a Pa-c Pc 


c^ + 2 



( 22 ) 


The conjugate momentum pT associated to matter appears linearly in H, and in this way 
a Schrodinger equation can be obtained with the matter variable playing the role of time. 
Therefore, all tools of ordinary quantum mechanics can, in principle, be employed in order to 
obtain predictions regarding the evolution of the universe. 

As usual, the Kantowski-Sachs model is quantized following Dirac’s quantization scheme 
by turning the variables pa, Pc and pT into operators, which fulfill the standard commutation 
relations 

. d ^ . d ^ . d 


Pa 


Pa = -T 


Pc 


Pc = 


PT 


PT = -1-. 


da dc dT 

In quantum cosmology, the Hamiltonian formulation of general relativity is employed through 
the ADM decomposition of the geometry and a Schrodinger-like equation, the Wheeler- 
DeWitt equation in the field representation, is constructed, which determines the wave function 
of the universe as a whole. Performing in this way, the Hamiltonian (22) becomes an operator. 
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which annihilates the physical states: H = 0 ^ = 0. Using the new time reparametrization 

t —>■ 2T we have the following Wheeler-DeWitt equation 


2 ^ , 2 _^ _ 2 
“ 9a2 + dc^ "" 



'I'(a, c,t) = 0 


(23) 


Functions in the classical phase space of this theory can commute. However, the associated 
operators defined above in Hilbert space, in general, do not commute. As a result, the order in 
which the product of two operators is written becomes significant, which leads to an ordering 
problem. The most direct way to write the operator is one that makes it symmetrical, which is 
a precondition for the operator to be self-adjoint, which will make it a physical observable. A 
short discussion about the self-adjointness of the operator is made in reference [2S]. Thus, we 
have 


We propose the following ordering 

o?dl -h qada = X . 

The inner product is, by definition, written as 

{ip, A(/>) = y V'* + qada^ (p da 

After some integrations by parts, we can write 


and similarly, if 
we will have 


(24) 

(25) 

(26) 


! Ip* {o?dl -h 2adc^ (p da = {Xip, cp) , 

(27) 

c^dl -)• c^dl -h qcdc = Y , 

(28) 

J Ip* (^c^d^ -|- 2 cdc^ (p dc = {Yip, cp) , 

(29) 


where qa = 2a, in equation (27), and qc = 2c, in equation (29), to ensure the hermiticity of H 
and the positivity of the Schrodinger inner product. A more detailed discussion can be found 
in [IB]. With this new ordering, the Wheeler-DeWitt equation can now be rewritten as 


da'^ 


d 

“da 


dc^ 


d 

dc 


a—1 


d 

dt 


'I'(a, c,t) = 0 


(30) 


Note that the signature of the kinetic term is hyperbolic, something that often occurs in 
anisotropic models. On the assumption of stationary solutions, 'k(a,c, t) = ip{a,c) we 

are led to 


2 9 ^ 2 2 < 9 ' 

C TTW + ^C- - c - a 

dc^ 


dc 


d 

da? ^ da 


i){a,c) = E 'ip{a,c) . (31) 


Analytical solutions of this equation are possible only for the particular case where a = 1 
(stiff matter). In addition to the reasons discussed above, this case is reasonable to check the 
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connection between Schutz’s formalism and cosmological models with anisotropic metrics. Thus, 
with a = 1 , we have 





(32) 


The operators + 2ada and + 2cdc in T^(0, oo) are unitarily equivalent to the operators 
9^ and d‘^ in L^(—oo,oo) m, which are self-adjoint (in fact, if two operators A and B are 
unitarily equivalent, B enjoy all the spectral properties that characterize A). This mapping is 
defined by the transformation 

y) = e“^) , 


where a = e ^ and 0 = 6^. Therefore, the equation (32) is modified to 

"^ 4 , y) = y) , 


(33) 


where we assume that E > {). 

Applying the method of separation of variables, we obtain 


that leads to the equations 


^{x,y) = X{x)Y{y) 
1 d'^X 


X dx‘^ 

where —is a separation constant and, 

1 d‘^Y 


= -A^ 


Y dy'^ 


- + e-^y) = 0 


(34) 

(35) 


(36) 


where = E — 

The general solution to the wave function of the universe follows from the above equations 


^Xa{x,y,t) = [Ai sm{Xx) + A 2 Cos{Xx)] [Bil„{e y) + B 2 K„{e ^)]e ^ ( 37 ) 

where Ai, A 2 , Bi and B 2 are integration constants, E and are modified Bessel functions 
and we assumed that a G C and A G M. 

Having this general solution, we need to investigate under what conditions the above wave 
function satisfies the boundary conditions, i.e., 'I'(a,c, t) = 0 (when a —>■ 0 and c —>■ 0). Under 
this restriction, we obtain the following wave function 

= . (38) 

With the boundary conditions, the wave packets may be constructed by taking superpositions 
of the A and a. The general structure of these superpositions is 

fOO fCO ,00 

^(x,y,t)=/ / C{X,a)e-^^^ ^dXda . (39) 

Jo Jo 





Note that the above solution admits the following factorization 


where 


Xa 


y, t) = Ia{x, t)Ix{y, t) , 

^0 

fOO ^ 

= / C2{a)K^{e-y)e-^^ Ua . 

Jo 


(40) 

(41) 

(42) 


For integral (41), if we adopt C'i(A) = e and X^t + Xx = {X^/t + to + “ 4 :{t+to) ’ 

we obtain 


A\/t~+~to 4" 


• /-OO 

( r 

X\ = e '‘(*+* 0 ) / 

cos 

Jo 

1 [ 


‘J\/t + to 


-.sm JjdA, 

(43) 

where to is an arbitrary constant. Introducing the variable u = X^/t^Jto + integral 

(43) takes the form 


Xa = 


g ®4(t+to) roo 


\/t + to 


J (cos — i sin du, 


(44) 


2 ^/t+t^ 


which has the solution given by 
. ^2 

g * 4 (t+tQ) ( 

Xa= , 1 - 2 C' 

\/t + to V 8 ( 

where C e S are the Fresnel integrals 


2 X 
TT 2A/t + to 


— z 


1-25 


2 X 

vr 2A/t + to 


(7(2;) = J cos ^9'^^ d9, 
S{z) = J sin^— 


For the probability density, we have 


|XaP = 

+ 5 


1 + 2 


4:(t + to) 

1 X 


c 


S 


1 X 
\/27r \/t + to 
1 X 


(K 


1 X 
\J2'K \/t + to 


-1 + 


(45) 

(46) 

(47) 

(48) 


\/27r -y/t + to / \ \ 'J2'K A/t + to 


- 1 


A similar procedure is employed to solve the integral (421. If we put the Bessel function in 
the integral form 

roo 

K^{c)= cosh iazdz, (49) 

Jo 

and we introduce (72(u) = and u± = a^t + to ± „ /m- ; obtain 


2 A/f++i- 

/tT (1 — z) —ccoshz+i ^-^ =- 

= e ^(V^)dz. 

V 8 ^/t + to JO 


(50) 
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Figure 1: The behavior of which is given by equation (48), and |Xo- 
equation (50). Here we adopt to = 1-5. 


which is found with 


The expressions found for \Tx\^ and |Xo-p are plotted in figure 

The relevant plots for the unnormalized probability density function is given in figures 
and|^ We note that the norm is time dependent and therefore the evolution is not unitary, which 
is also found in the context of the Bianchi model [5U] . While this result may reinforce evidence of 
the existence of an incompatibility between anisotropic models and Schutz formalism, it should 
be noted that, according to the equation (33) the operator associated with the wave function 
obtained is self-adjoint. Since the initial conditions are also well established, the loss of unitarity 
must have been caused by the wave packet chosen. Although the packet is nothing more than 
a superposition of plane wave functions, it is possible that certain choices are not able to locate 
the solution in a finite region of space or to preserve the norm independent of time. This may 
explain the occurrence of pathologies observed in this sense in quantum cosmological models. 
The wave packet could have been built someway else, which in this case means other choices 
for C'i(A) and 6*2 (<t). The choices made here allowed us to find an analytical expression for the 
wave function that although useful, may not lead to the best physical description. In our case 
the price to pay is a unnormalizable wave packet, since the normalization integral diverges in 
the X direction. 


3 Conclusions and Discussions 

In this paper, we studied the Kantowski-Sachs quantum cosmological model. The introduction 
of a time variable in the theory was performed using the Schutz’s formalism, which associates 
the time with degrees of freedom of a perfect fluid coupled to gravitation, chosen in this work 
as stiff matter. 

The Schutz’s formalism enables the introduction of time in theory to occur consistently, al¬ 
lowing the Wheeler-DeWitt equation in this minisuperspace to have an explicitly time-dependent 
solution. However, the model is not unitary since the norm is time-dependent, which has been 
observed in the context of Bianchi I model [20]. The guarantee that the operator associated 
with the model is self-adjoint was not able to eliminate anomalies already detected in anisotropic 
cosmological models such as the loss of unitarity. This indicates that certain conditions may 
arise at later stages of the model construction, as in the build procedure of the wave packet. In 
fact, even if we have built a package with conventional auxiliary functions (such as Gaussian 
type functions), it is possible that it deeply affects the behavior of certain solutions, leading to 
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Figure 3: The wave packet as a function of a and t, for c = 0, 5, 10 and to = 1.5. 


the problems encountered in the context of this model and others. In this sense, although Schutz 
formalism is not necessarily responsible for anomalous aspects of the operator that describes the 
quantum system, its use may cooperate to make the subject solution behavior changes during 
the procedure of obtaining the package. 
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Figure 4: The wave packet as a function of a and c, for t = 0, 5, 10 and to = 1-5- 

The possibility that certain functions used in the construction of the wave packet can lead 
to the appearance of anomalies and physically inadequate results should be investigated more 
thoroughly in future work. An alternative is to compare solutions based on different wave 
packets. Other perspectives are the construction of models with other types of fluid or even 
the introduction of a second fluid (though this is a somewhat artificial process), which would 
further increase the number of degrees of freedom of the minisuperspace. This may allow a 
deeper understanding of the connection between matter and geometry in the early universe. 
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